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Abstract 

Background: More than two-fifths of the world's population uses solid fuels, mostly biomass, for cooking. The resulting 
biomass smoke exposure is a major cause of chronic obstructive pulmonary disease (COPD) among women in developing 
countries. 

Objective: To assess whether lower woodsmoke exposure from use of a stove with a chimney, compared to open fires, is 
associated with lower markers of airway inflammation in young women. 

Design:\Ne carried out a cross-sectional analysis on a sub-cohort of participants enrolled in a randomized controlled trial in 
rural Guatemala, RESPIRE. 

Participants: \Ne recruited 45 indigenous women at the end of the 18-month trial; 19 women who had been using the 
chimney stove for 18-24 months and 26 women still using open fires. 

Measurements:VJe obtained spirometry and induced sputum for cell counts, gene expression of IL-8, TNF-a, l\/ll\/lP-9 and 12, 
and protein concentrations of IL-8, myeloperoxidase and fibronectin. Exhaled carbon monoxide (CO) and 48-hr personal CO 
tubes were measured to assess smoke exposure. 

Resuits:MMP-9 gene expression was significantly lower in women using chimney stoves. Higher exhaled CO concentrations 
were significantly associated with higher gene expression of IL-8, TNF-a, and l\/ll\/lP-9. Higher 48-hr personal CO 
concentrations were associated with higher gene expression of IL-8, TNF- a, l\/ll\/lP-9 and IVlMP-12; reaching statistical 
significance for MMP-9 and MMP-M. 

Conclusions: Compared to using an open wood fire for cooking, use of a chimney stove was associated with lower gene 
expression of l\/ll\/lP-9, a potential mediator of airway remodeling. Among all participants, indoor biomass smoke exposure 
was associated with higher gene expression of multiple mediators of airway inflammation and remodeling; these 
mechanisms may explain some of the observed association between prolonged biomass smoke exposure and COPD. 
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Introduction 

Household air pollution from cooking with solid fuel (wood, 
dung, crop residues, coal, charcoal) is estimated to be responsible 
for 4.6% of worldwide disability adjusted life-years lost. [1] This is 
due to acute lower respiratory infections (ALRI) in children, 
chronic obstructive pulmonary disease (COPD), lung cancer, 
cataracts, and cardiovascular disease in adults. [2] More than two- 
fifths of the world's population uses solid fuel as their primary fuel 
for cooking and the vast majority of these people do so over open 



fires indoors. [3] This inefficient way of cooking exposes women 
(who tend to perform the majority of cooking) on a daily basis to 
high levels of smoke. Mean indoor 24-hour concentrations of fine 
particulate matter (particles with aerodynamic diameter of <2.5 
microns, PM2.3), a component of biomass smoke, have been found 
to be up to 100 times higher than the air quality guidelines of the 
World Health Organization (10 ug/m3 annual average) and are 
commonly 30 times higher [4,5]. 

Epidemiologic studies of women from developing countries 
have shown strong associations between use of biomass cooking 
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fuel and chronic bronchitis and/ or COPD. [6-15] Separate meta- 
analyses by Hu, et al. and Kurmi, at al identified an increased risk 
of COPD among those exposed to biomass smoke with odds ratios 
of 2.44 (95% CI, 1.9-3.33) and 2.80 (95% CI, 1.85-4.0), 
respectively. [16,17] A prospective study showed that participants 
with COPD due to biomass smoke have similar mortality rates as 
participants with COPD due to tobacco smoke [18]. 

Multiple studies have identified an increase in inflammatory 
mediators in sputum that may suggest a mechanism for this 
increased risk. One study of COPD due to biomass smoke 
exposure sampled respiratory tract lining fluid by bronchoalveolar 
lavage (BAL), finding increased matrix metaUoproteinase (MMP)- 
12 gene expression and activity and increased MMP-9 activity. 
[19] More recently, observational studies of individuals in West 
Bengal, India exposed to a combination of solid fuel types (wood, 
dung, and crop residues) found an increase in inflammatory cells 
and other markers of inflammation in expectorated sputum, 
including total cell counts, leukocyte counts, and sputum levels of 
interleukin (IL)-6, IL-8, tumor necrosis factor (TNF)-0(, and 
reactive oxygen species after adjustment for confounders, though 
this study was limited by significant differences between the 
biomass smoke-exposed and comparison (liquid petroleum gas 
stove) groups at baseline. [20-22] Our study has several 
advantages over previous literature on this topic, as subjects were 
recruited following a randomized stove intervention trial and more 
intensive exposure monitoring was collected on participants 
including exhaled breath CO and 48-hour personal CO concen- 
trations. 

We carried out our study in rural highland Guatemala, the site 
of a randomized controlled trial (Randomised Exposure Study of 
Pollution Indoors and Respiratory Effects - RESPIRE) in which 
households using open fires indoors for cooking were randomly 
assigned to receive a well-operating chimney stove (locally called 
the Plancha), either at the start or end of the study, after an 18- 
month study period. [23] The main outcomes of RESPIRE were 
incidence of acute lower respiratory infections in infants and 
respiratory symptoms and lung function in their mothers. 
Measurements showed that the chimnc^y stoves, which were 
indigenously developed in Guatemala and have been deployed 
widely since the late 1980s, were popular, sturdy, and redu[:ed 
kitchen air pollution and personal exposures to women and 
chfldren [24,25]. 

We used induced sputum to obtain samples of respiratory tract 
lining fluid from a subgroup of young women enrolled in 
RESPIRE at the end of the 18-month trial, to determine if 
reduced smoke exposure from use of the chimney stove is 
associated with decreases in markers of airway inflammation and 
injury that characterize COPD. [26-30] We chose to examine 
biomarkers suggestive of macrophage and epithelial cell activation 
(IL-8, TNF-a, MMP-9 and 12), neutrophilic inflammation 
(differential cell counts, myeloperoxidase and MMP-9), and injury 
(fibronectin), each of which have been linked to COPD 
pathogenesis in previous studies. 

Methods 

Ethics Review 

The study protocol was approved by the institutional review 
boards at the University of California, San Francisco; the 
Universidad del VaUe, Guatemala; the University of Bergen, 
Norway; and the University of Liverpool, United Kingdom. 
Because of the low literacy rate in this population, oral informed 
consent was obtained in each participant's native language and 
documented by field workers using a standardized form. This 



method of consent was approved by each institutional review 
board. 

Study Design and Participants 

RESPIRE was a randomized controlled stove intervention trial 
undertaken between October 2002 and December 2004 in which 
households with children up to 18 mo of age were randomized 
either to receive a stove with a chimney or to continue cooking 
over an open fire until the end of the trial. Details of the RESPIRE 
protocol have been published elsewhere. [31,32] Participants of 
this cross-sectional study were all mothers of children participating 
in RESPIRE, recruited in December 2004 near the end of the 18- 
month study period. A convenience sample of 50 women was 
invited to participate in this cross-sectional study and 45 accepted, 
with 26 from open-fire and 19 from chimney-stove homes. AH tests 
were done at the Ministry of Health clinic in San Lorenzo over a 
2 -week period, where participants underwent spirometry and 
sputum induction to assess for airway markers of inflammation 
and injury. 

Women were assessed by self-report for the following exclusion 
criteria: a) acute respiratory infection within the past 4 weeks; b) 
use of steroids; c) history of current or former smoking; or d) lack of 
participation in cooking. One (2%) woman was excluded because 
of acute respiratory infection. 

Spirometry 

Using a portable spirometer (EasyOne, ndd), we measured 
forced expiratory volume in 1 second (FEVi) and forced vital 
capacity (FVC) before and 15 minutes after administration of a 
bronchodilator. [33] Only measurements that were acceptable and 
reproducible by 1994 American Thoracic Society criteria were 
analyzed. [34,35] Reversibifity was defined as S12% and S 
200 ml increase in FEVi or FVC, and COPD as post-broncho- 
dilator FEVi/F\^C ratio <0.70 [36]. Reference values were 
calculated from the NHANES Ill-derived equation for Mexican- 
Americans. [34] Concavity of expiratory flow-volume curves was 
assessed by a puknonologist blinded to the RESPIRE stove 
intervention status. 

Sputum Induction and Processing 

One of the investigators induced and processed the sputum and 
performed total cell counts at the Ministry of Health clinic with the 
assistance of a RESPIRE field worker. The cell slides and frozen 
sputum samples (supernatant and cells) were transported to the 
Lung Biology Center at the University of California, San 
Francisco. 

Cell Counts, Gene Expression and Protein Concentrations 

Differential cell counts (400 non-squamous cells) were per- 
formed in duplicate by two independent counters who were 
blinded to the RESPIRE stove intervention status. All gene 
expression and protein assays were conducted at the Lung Biology 
Center by one of the investigators. Concentrations of myeloper- 
oxidase (MPO), IL-8 and fibronectin in sputum supernatants were 
measured using commercially available ELISAs (MPO: Calbio- 
chem, Germany; IL-8: R&D Systems, Minneapolis, USA; 
fibronectin: Biomedical Technologies, Massachusetts, USA). 
Relative quantification of IL-8, TNF-a, MMP-9 and MMP-12 
gene expression was performed by comparison to an endogenous 
control, glyceraldehyde 3-phosphate dehydrogenase (GAPDH), 
using the two-step Real-Time Quantitative Reverse Transcription- 
Polymerase Chain Reaction (Applied Biosystems, TaqMan) 
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according to standard protocols available at http://docs. 
appliedbiosystems.com/pebiodocs/04371095.pdf). 

Personal 48-hr Carbon Monoxide (CO) 

Per.sonal 48-hr CO concentrations [reported in parts per million 
(ppm)] were measured using Drager Safety (Luebeck, Germany) 
passive difiFusion tubes, which were carried by participants in small 
loosely woven bags previously shown not to impede difiFusion of 
CO. [24] The passive diffusion tubes function by changing colors 
as CO passes through the tubes, with the length of the color stain 
being proportional to total CO exposure in ppm over the 48-hour 
period. Exhaled breath CO and personal 48-hr CO monitoring 
were used in this study as surrogates for biomass smoke exposure; 
Personal CO been shown to be a reliable surrogate for PM2.5 
exposure due to biomass combustion. [37] Quality control and 
calibration of personal CO measurements were maintained 
according to previously described protocols. [4] These were 
collected at least two times after stove installation (generally 0, 3, 6, 
and 12 months). As the number of personal CO measurements 
and timing relative to sputum induction was variable, the last 48- 
hr personal CO measurements obtained prior to sputum induction 
were used for regression analyses. 

Exhaled Carbon Monoxide Measurement 

During the 18-month duration of the RESPIRE study, exhaled 

CO concentration (in ppm) was measured for each participant 
using the Micro CO Meter (Micro Medical) at three separate times 
(6, 12 and 18 months). Exhaled CO measurements were not 
obtained during the induced sputum visits. The maneuver 
required the participant to hold her breath for 15 seconds and 
exhale completely with her lips sealed around a disposable 
mouthpiece. At each testing session, the participant performed 
the maneuver three times and the two highest readings were used 
in analysis. The mean value of the six highest readings obtained 
during the three testing sessions was used as the predictor variable 
in exposure-responses analyses. Two women were unable to 
complete any exhaled CO measurements. Four women completed 
two of three exhaled CO measurements; their mean CO was 
calculated as the mean of the two highest reading taken at each 
measurement. Exhaled CO monitors were calibrated against a 
standard CO span gas every 2-4 weeks during the study period. 

Statistical Analysis 

Statistical analysis was performed using Stata software version 
12.1 and SAS software version 9.1.3. A sample size of 17 per 
group was calculated to detect a difference of 10% in percent 
neutrophils between the groups, with type 1 error of 0.05 (two- 
sided) and power of 80%. Univariate comparisons between the 
two groups were performed using unpaired student t-Tests for 
normally distributed continuous variables or their logarithmic 
transformations, WUcoxon rank sum for non-normally distributed 
continuous variables, and chi-square or Fisher exact tests for 
categorical variables. Normally distributed data are expressed as 
mean values with standard deviation (SD) and non-normally 
distributed data are expressed as median values with interquartile 
ranges (IQR). 

Simple linear regression models were used to describe exposure- 
outcome associations, with untransformed exhaled CO concen- 
tration (ppm) and 48-hr personal CO concentrations as the 
independent variable and the inflammatory markers or their 
logarithmic transformations as dependent variables; Simple linear 
regression was also used to describe the relationship between lung 
function and percent sputum neutrophils. Distributions were 
graphically assessed and variables were log-transformed for 



regression analysis when non-normally distributed (resulting in 
log transformation of all gene products). Assumptions for linear 
regression were assessed using standardized residual versus fitted 
value and normal rjuantile plots (to assess functional form and 
constant variance assumptions) and Shapiro WHk tests were used 
to formally assess for normahty of residual distributions; assump- 
tions were met for all included regression analyses. 

Results 

Participant Characteristics 

The characteristics of women, at the time of induced sputum 
visits, from open-fire and chimney-stove homes are shown in 
Tables 1 and 2. There were no significant differences in age, 
height, weight, village, highest educational attainment (for 
participant or spouse), primary cooking location, or in the 
following pre-or post-bronchodilator spirometric indices: percent 
predicted FEVi, percent predicted FVC, and FEVi/FVC ratio. 
Although the percentages of women with no formal education and 
with their primary cooking area in the same room as the bedroom 
were higher in the open fire group, these difierences were not 
significant. Thirty-five (78%) women had acceptable (as deter- 
mined by American Thoracic Society performance criteria) pre- 
and post-bronchodilator spirometric measurements, an additional 
six (13%) women had acceptable pre- or post-bronchodilator 
spirometric measurements, and four (9%) women had unaccept- 
able pre- and post- measurements. Most (16 out of 19 women, or 
84"'o) had the chimney stove for at least 18 months, with the 
remaind(-r having their stove at least 15 months. The participants 
from the open-fire group did have more cases of both concavity of 
expiratory flow volume curve (50% versus 30%, p = 0.22) and 
airflow reversibility (16% versus 0%, p = 0.23), although these 
differences were not statisticaUy significant. Though no participant 
was known to have respiratory disease at the time of enroUment, 
one participant was identified as having COPD based on 
spirometry performed during the study. 

Exposure Characteristics 

Exhaled CO measurements were available for 43 women and 
48-hr personal CO measurements were available for all partici- 
pants. Participants from the open-fire homes (n = 26) had 
significantiy higher mean exhaled CO concentrations than women 
from chimney-stove homes (n= 17), with 7.1 ppm (SD 1.6) and 
5.2 ppm (SD 1.1), respectively, with a mean difference of 
1.96 ppm (95% CI 1.1-2.8, p = 0.0003). Women from open fire 
homes (n = 26) also had significantly higher 48-hr personal CO 
concentrations than women from chimney-stove homes (n = 19), 
3.5 ppm (SD 1.7) and 1.5 ppm (SD 1.2), respectively, with a mean 
difference of 2.0 ppm (95% CI 1.0-2.9, p<0.001). Pearson's 
correlation coefficient for exhaled CO and 48-hr personal CO was 
0.48. 

Induced Sputum Quality and Cell Counts 

Of the 45 women eligible for sputum induction, 40 (89%) 
women produced sputum and 31 (69%) women had cell slides of 
acceptable quality for use in the analysis: 14 from open-fire and 17 
from chimney-stove homes. Total and differc-iitial cell counts are 
shown in Table 3. There were no significant difierences between 
the groups in total leukocytes, differential counts of neutrophils, 
macrophages, eosinophils and lymphocytes. Percent epithelial ceUs 
was significantiy higher in participants from the open-fire group in 
comparison to the chimney-stove group (p = 0.05). In post hoc 
analysis, percent neutrophils in sputum was negatively associated 
with percent predicted FEVi (P = -0.26, p = 0.01; Figure 1). 
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Table 2. Lung function of study participants. 





Pre-bronchodilator spirometric indices 


Cliimney stove n = 1 7 


Open fire n = 20 


p-value 


95%CI for difference 


FEVl, percent predicted (mean, SD) 


102 (13) 


103 (12) 


0.76 


-7, 10 


FVC, percent predicted (mean, SD) 


108 (13) 


111 (13) 


0.4 


-5, 12 


FEVl/FVC % (mean, SD) 


0.82 (0.07) 


0.80 (0.05) 


0.54 


-0.06, 0.03 


"Concavity, n (%) 


5 (30) 


10(50) 


0.22 




Post-bronchodilator spirometric diagnostics 


n = 16 


n = 19 






"Reversibility, n (%) 


0 


3 (16) 


0.23 




•COPD, n {%) 


1(6) 


0 


0.46 




Pre-bronchodilator spirometric indices 


CKiimney stove n = 1 7 


Open fire n = 20 


p-value 


95%CI for difference 


FEVl, percent predicted (mean, SD) 


102 (13) 


103 (12) 


0.76 


-7, 10 


FVC, percent predicted (mean, SD) 


108 (13) 


111 (13) 


0.4 


-5, 12 


FEVl/FVC % (mean, SD) 


0.82 (0.07) 


0.80 (0.05) 


0.54 


-0.06, 0.03 


•Concavity, n (%) 


5 (30) 


10(50) 


0.22 




Post-broncliodilator spirometric diagnostics 


n = 16 


n = 19 






"Reversibility, n (%) 


0 


3 (16) 


0.23 




*COPD, n (%) 


1(6) 


0 


0.46 





*Data are categorical variables and the percent represents number of participants with concavity or reversibility divided by total number of participants in each group. 
doi:1 0.1 371 /journal.pone.0088455.t002 



Induced Sputum Protein Concentrations 

Of the 40 women who produced sputum, 32 (80%) women 
produced enough sputum for protein analysis. Of these, 22 (69%) 
women had samples that were sufficient in quantity to carry out all 
three of the protein assays. Mean protein concentrations of IL-8 
were lower in the chimney-stove group (1801 ng/ml versus 
2573 ng/ml, p = 0.15) while mean fibronectin (485 ng/ml versus 
355 ng/ml, p = 0.44) and median MPO (78 ng/ml versus 26 ng/ 



ml, p = 0.40) concentrations were higher in this group, though 
none of these difiFerences reached statistical significance (Table 4). 

Induced Sputum Gene Expression 

Good quality mRNA is difficult to extract from induced sputum 
because samples tend to have few cells and these difficulties may 
be enhanced when initial processing is carried out at a rural field 
station. [38] Despite these challenges, we recovered RNA from die 
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Figure 1. Simple linear regression between percent neutrophils (independent variable) and percent predicted FEV^ (dependent 
variable). Dots represent individual estimates. Solid line is made to linear best fit model. 
doi:1 0.1 371 /journal.pone.0088455.g001 
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Table 3. Induced sputum cell counts of study participants. 





Cell concentrations 


Chimney stove 


Open fire 


p-value 




(n = 14) 


(n = 17) 




Total leukocytes x 10'*/ml {mean, SD) 


234.4 (171) 


255.2 (172) 


0.7 


Neutrophils x 10*/ml (mean, SD) 


187.8 (167) 


189.6 (174) 


0.71 


Macrophages x lO'^/ml (mean, SD) 


54.9 (52) 


88.7 (94) 


0.61 


Differential cell counts 


^ Neutrophils, % WBC (median, IQR) 


79.3 (59.4, 87.3) 


72.5 (48.5, 90.6) 


0.58 


^Macrophages, % WBC (median, IQR) 


20.7 (12.7, 40.1) 


27.5 (9.2, 51.5) 


0.63 


*Epithelial cells, % all (median, IQR) 


1 .25 (0.30, 2.0) 


2.0 (1.3, 3.3) 


0.05 



^The mean of cell type divided by the sum of all white blood cells multiplied by 100. 
■'"The mean of epithelial cells divided by sum of all cells multiplied by 100. 
doi:l 0.1 371 /journal.pone.0088455.t003 



34 samples available for gene expression analysis and 29 samples 
(85%) were of good quality. Three (9%) assays were rejected 
because duplicate Ct measurements differed by more than 0.5. 
MMP-12 and -9 gene expression was low in the samples - all 
samples had mean Ct levels above 28. As the accuracy of gene 
expression measurements at such low concentrations can be 
limited, we carried out our analyses both including and excluding 
samples with very low levels of gene expression (defined as mean 
Ct levels >37) and found similar results (n = 22 for MMP-9, n = 27 
for MMP-12). We used the Comparative Ct Method to quantify 
our gene of interest as the number of gene copies per 10* copies of 
GAPDH because validation experiments revealed efficiencies of 
GAPDH, IL-8, TNF-alpha and MMP-9 to be above 90%. The 
chimney-stove group had significantly less MMP-9 gene expres- 
sion than the open-fire group (0.002 vs 0.01, p = 0.01). In addition, 
gene expression of TNF-alpha was more than 50% lower and gene 
expression of MMP- 1 2 was 1 0 fold lower in the chimney stove 
group, though these differences were not significant (Table 5). 

Exposure-Outcome Associations 

Exhaled CO concentrations (in ppm) were positively associated 
with MMP-9 (B = 2.17, 95% CI 1.2-3.9), TNF-a (B = 1.3, 95% CI 
1.0-1.7) and IL-8 (B= 1.64, 95% CI 1.1-2.4) gene expression in 
log copies per 10* copies of GAPDH (Table 6). Personal 48-hr CO 
concentrations (in ppm) were positively associated with MMP-9 
(B = 2.58, 95% CI 1.5-4.4) gene expression. Additionally, there 
were positive, non-significant associations between 48-hr personal 
CO exposure with MMP-12 (B = 0.4, 95%CI 0.96-2.0), IL-8 



(B= 1.41, 95% CI 0.98-2.0) and TNF-a (B= 1.22, 95% CI 0.98- 
1.5) gene expression. There were no significant relationships 
between concentrations of either exhaled CO or personal 48-hr 
CO exposure and fibronectin, IL-8, MPO protein concentrations, 
percent neutrophils, or continuous spirometric variables (% 
predicted FEVl, % predicted FVC, or FEVl/FVC ratio). 

Discussion 

We took advantage of a unique opportunity - a randomized 
stove intervention trial - to compare inflammatory markers in 
sputum between difFerentially woodsmoke exposed, young, indig- 
enous Guatemalan women. Our cross-sectional analysis shows that 
women from the chimney-stove homes had significantly fewer 
epithelial cells and lower MMP-9 gene expression in induced 
sputum when compared to women from open-fire homes. 
Exposure-outcome analyses demonstrated significantly higher 
levels of IL-8, TNF-ot, MMP-12 and MMP-9 gene expression 
associated with higher levels of biomass smoke exposure, using 
personal CO and exhaled breath CO as surrogates for exposure. 
Additionally, the direction of the association between all examined 
gene expression products and these surrogates for biomass smoke 
exposure was positive, with higher expression of all gene products 
at higher levels of biomass exposure, with the exception of 
MMP 12 and exhaled CO. As the protein products of these genes 
are associated with COPD pathogenesis in tobacco smokers, these 
findings suggest that important commonalities may exist between 
tobacco and biomass smoke-associated COPD pathogenesis. 



Table 4. Induced sputum protein concentrations of study participants. 





Protein concentration 


Chimney stove 


Open fire 


p-value 


IL-8, ng/ml (mean, SD) (n = 32) 


1801 (1406) (n=15) 


2573 (1533) (n= 17) 


0.15 


Fibronectin, ng/ml (mean, SD) (n = 22) 


485 (409) (n=11) 


355 (374) (n=n) 


0.44 


Myeloperoxidase, ng/ml (median, IQR) (n = 29) 


78 (16, 625) (n= 15) 


26 (14, 463) (n= 14) 


0.4 


Protein concentration 


Chimney stove 


Open fire 


p-value 


IL-8, ng/ml (mean, SD) (n = 32) 


1801 (1406) (n= 15) 


2573 (1533) (n= 17) 


0.15 


Fibronectin, ng/ml (mean, SD) (n = 22) 


485 (409) (n=11) 


355 (374) (n=n) 


0.44 


Myeloperoxidase, ng/ml (median, IQR) (n = 29) 


78 (16, 625) (n= 15) 


26 (14, 463) (n=14) 


0.4 



doi:1 0.1 371 /journal.pone.0088455.t004 
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Table 5. Induced sputum gene expression of study participants. 



Gene copies per 1 o" copies of GAPDH ^ 


Chimney stove 


Open fire 


p-value 


TNF-a (median, IQR) 


299 (111, 509) 


682 (165, 1566) 


0.12 


(n = 29) 


(n = 14) 


(n=15) 




IL-8 (median, IQR) 


519 (20, 2016) 


575 (268, 1 349) 


0.55 


(n = 29) 


(n = 14) 


(n=15) 




MMP-12 (median, IQR) 


0.01 (0.005, 0.08) 


0.10 (0.02, 0.22) 


0.17 


(n = 27) 


(n = 14) 


(n = 13) 




MMP-9 (median, IQR) 


0.002 (0.0008, 0.006) 


0.01 (0.003, 2.96) 


0.01 


(n = 22) 


(n=12) 


(n=10) 





^glyceraldehyde 3-phosphate dehydrogenase. 
doi:l 0.1 371 /journal.pone.0088455.t005 



particularly if these findings are replicated and extended in larger 
longitudinal studies. 

Those with COPD due to tobacco smoke have been shown to 
have increased levels and activity of matrix metalloproteases 
(MMP-8, -9, and 12) as well as IL-8 and TNF-ot in induced 
sputum and BAL fluid. [39-47] IL-8 and TNF-a are potent 
cytokines released by activated epithelial cells and macrophages. 
MMPs are proteases that are involved in tissue growth, injury and 
repair as well as modulating inflammatory cascades. [48] To our 
knowledge, only one previous study has examined MMP levels in 
women with COPD due to biomass smoke and found increased 
levels of MMP-9 activity and increased MMP— 12 gene expression 
and activity in BAL when compared to healthy controls. [19] 
Although sputum induction samples both proximal and distal 
airway lining fluid, our data are consistent with and extend these 
fmdings. While most of our participants had normal spirometry, 
the prevalence of mild airflow obstruction (concave expiratory 
flow-volume curves) in these young women does suggest the 
potential for the development of woodsmoke-induced COPD. 
Given the observed associations between smoke exposure and 
gene expression identified in this study, reductions in smoke 
exposure with a cleaner stove and/or fuel type may prevent the 
chronic pulmonary inflammation and airway remodeling that 



culminates in COPD in many such biomass smoke exposed 
women. 

This study has several noteworthy limitations. Given the cross- 
sectional nature, we cannot say with certainty that the observed 
differences in gene expression are due to the chimney stove, as 
baseline levels for these biomarkers were not collected and the 
identified abnormahties could have preceded the intervention. 
There is some risk of selection bias since we employed convenience 
sampling; while participants were largely selected based on 
proximity to the field office, it is possible that participants with a 
larger reduction in indoor smoke, symptoms, or both may have 
chosen to participate. This sampling strategy also limits the 
generalizability of these results as this sample may not adequately 
represent biomass smoke-exposed women both inside and outside 
of Guatemala. As exposures were not collected immediately prior 
to outcome measures and only represent short periods of exposure, 
there is a risk of non-differential exposure misclassification, which 
most frequendy results in an underestimate of the true effect size. 
Additionally, while gene expression was elevated for multiple 
genes implicated in COPD pathogenesis, the protein products or 
activity of several of these genes were either not measured (for 
MMP-9, MMP-12, and TNF-ot) or were not convincingly 
associated with biomass smoke exposure (for IL-8). Concordant 
results for gene expression and protein product analyses would 



Table 6. Exposure-outcome relationships between exhaled concentration of CO and inflammatory markers. 



Gene copies per 10" Exhaled Breath CO (per 1 ppm) Personal CO Tube (per 1 ppm) 



copies of GAPDH^ 


Coefficient 


95% CI (p-value) 


Coefficient 


95% CI (p-value) 


MMP-9 (n = 22) 


2.17* 


1.2, 3.9 (0.01) 


2.58* 


1.5, 4.4 (0.01) 


TNF-a (n = 29) 


1.30* 


1 .0, 1 .7 (0.04) 


1.22* 


0.98, 1 .5 (0.08) 


IL- 8 (n = 29) 


1.64* 


1.1, 2.4 (0.02) 


1.41* 


0.98, 2.0 (0.06) 


MMP-12 (n = 27) 


0.93* 


0.6, 1 .3 (0.84) 


1.49* 


1.0, 2.2 (0.05) 


Protein concentrations ng/ml 


Fibronectin 


-8.6 


- 1 28, 1 1 1 (0.88) 


-63.5 


-152, 25.3 (0.15) 


Myeloperoxidase 


4.3 


-88.6, 97.3 (0.92) 


-26.9 


-104, 50.7 (0.48) 


IL-8 


170 


- 1 34, 474 (0.26) 


49.3 


-245, 343 (0.73) 


Differential cell counts 


Neutrophils, %WBC 


0.2 


-5.1, 5.5 (0.93) 


1.11 


-4.5, 6.6 (0.68) 



'glyceraldehyde 3-phosphate dehydrogenase. 
*Log-linear scale after back-transformation. 
doi:l 0.1 371/journal.pone.0088455.t006 
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strengthen the observations in this study. Finally, our assessment of 
multiple outcomes increases the risk of a Type I error, though the 
consistency of the findings serves to diminish this concern to some 
degree. 

Unlike recent observational studies in West Bengal, which 
reported lower numbers of airway neutrophils among individuals 
using LPG stoves compared to biomass stoves, we did not identify 
a similar finding in the chimney-stove group in our study. [20,22] 
The persistence of high airway neutrophils in the group receiving a 
chimney stove may be explained by the long historical exposure to 
woodsmoke of our participants, rather than their more recent 
exposure. Assuming relatively stable time-activity patterns in this 
young, homogenous group of women, we examined the associa- 
tion between age and sputum percent neutrophils. Although the 
mean age of our participants was only 26 years, our analysis did 
show a trend towards an association between age and percent 
neutrophils (P=1.4 p = ().()7. Data not shown). If cumulative 
exposure determines nc'utrophilic inflammation, it may take more 
time (i.e., more than the 18-month time period of RESPIRE) with 
reduced exposure to see a reduction in airway neutrophils. 

Alternatively, as there was a notable overlap in exposure profiles 
in the chimney-stove and control groups despite significant 
reductions in mean personal exposures, it is possible that a 
threshold effect exists for airway inflammation due to biomass 
smoke that was exceeded by the majority of individuals in both 
groups. We know that women from the chimney-stove group 
continue to have woodsmoke PM exposure much higher than 
recommended by WHO Air Quality Guidelines. Collectively, 
these findings could have important implications for future stove 
intervention trials focused on reducing COPD incidence or 
progression due to biomass smoke exposure, suggesting that 
longer periods of follow-up and cleaner stove interventions may be 
necessary to more clearly demonstrate an effect. 

Exhaled CO measurements have been used repeatedly to 
distinguish between smokers and non-smokers (cut-ofFs between 6- 
8 ppm) and recently, in a study from Taiwan, Hung et al. showed 
an association between increasing tobacco consumption and 
exhaled CO concentration. [49] We used exhaled CO to assess 
personal exposure to woodsmoke because for cultural and 
economic reasons tobacco smoking is rare in both men and 
women in our study area. We used average values from three time 
points as a measure of chronic exposure, although the value from 
each time point is reflective of acute exposure. The addition of 
personal 48-hour CO measurements, which should capture 
exposure from multiple stove burn cycles and is strongly correlated 
with PM2.5 exposure in this population, further strengthens the 
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